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A new design for heavy-ion beam driven isentropic compression experiments is suggested and 
analysed. The proposed setup utilises the long stopping ranges and the variable focal spot geometry 
of the high-energy uranium beams delivered at the GSI Helmholtzzentrum fur Schwerionenforschung 
and Facility for Antiproton and Ion Research accelerator centers in Darmstadt, Germany, to pro- 
duce a planar ramp loading of various samples. In such experiments, the predicted high pressure 
amplitudes (up to 10 Mbar) and short timescales of compression (below 10 ns) will allow to test the 
time dependent material deformation phenomena at unprecedented extreme conditions. 

PACS numbers: 62.50.-p, 64.60.A-, 52.50.Gj 



Research carried out in astrophysics, planetary and 
material sciences seek a thorough understanding of the 
behavior of matter at high pressures. For example, the 
equation of state (EOS) of iron under pressures of 1- 
4 Mbar is crucial in order to determine the state of the 
Earths core [JJ 03 [3]. EOS data around 0.1 Mbar is 
required to establish the state and the composition in 
Earths lower mantle |3J H], while the dynamics of the 
processes in the mantle is considered to be dependent 
on the structural phase transformation kinetics [5j. The 
pursuit after materials for technological applications [B] 
also entails a detailed understanding of the kinetics of the 
high- pressure phase transitions. Modelling the physical 
processes during the projectile impact relevant to mc- 
tcoroid protection and crater formation [7] requires the 
dynamical response of solids at ultrahigh strain rates. 
Hydrogen EOS at high pressures are especially impor- 
tant for understanding the structure and evolution of the 
hydrogen-bearing astrophysical objects such as the giant 
planets like Saturn and Jupiter [8j [9] as well as for the 
inertial fusion energy research [TUJ [TT] . 

The technique of isentropic compression using ramp 
wave loading (RWL) [Tl[Tl[Tl[Tl[Tl[T71[THl[ini[2ni 
|2"T1 |2"21 |2"5"] allows one to sample EOS data along an isen- 
trope up to several Mbar pressures, located on the phase 
diagram between the parameter regions accessible in di- 
amond anvil cell 24, 25j and shock wave (26j [27] exper- 
iments. Unlike a shock wave experiment where a single 
point on a shock adiabat is obtained, in isentropic com- 
pression experiments a continuous set of data points is 
recorded and the solid state of a sample is ensured up 
to high pressures. RWL technique was also shown to be 
a more sensitive tool for studying the dynamics of ul- 
trafast structural phase transformations than the shock- 
wave based [251 HI [30] methods. 

RWL has been demonstrated with different drivers, 
such as magnetic pulse loading using high-current pulsed 
power generators O [13j [14j [15j [TBI Hl]j high-power 
lasers [20j EU E21 [23] , and gas guns [18] or high explo- 
sives [JJ5] using graded density impactors. The typical 



loading times with these drivers are 10 ns, 100 ns and 
1 [is, respectively. In this Letter, a new scheme for pla- 
nar isentropic compression experiments using an intense 
heavy ion beam as a driver is proposed and analysed. 
The long absorption range and the variable focal spot 
size and shape of energetic ion beams allows one to de- 
sign isentropic compression experiments with fairly pla- 
nar geometry. The beam parameters needed to generate 
pressures of up to 4 Mbar in aluminium and 7 Mbar in 
iron are well within the reach of the uranium beams that 
will be delivered at the new Facility for Antiproton and 
Ion Research (FAIR) which is being built in Darmstadt. 
The beam intensities needed to generate pressures ap- 
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FIG. 1: Schematic of the experimental design. The ion beam 
with an elliptic focal spot propagates in z direction. It heats 
the absorber slab along its surface. The depth of the beam 
inside the absorber, i.e. the distance between the beam's axis 
and the inner surface of the slab is Ar. The target foil of 
thickness d is placed parallel to the absorber at a distance D. 
Hatched areas represent the support washer used to and to 
coax the expanding absorber material. The velocity of the 
rear target surface is recorded using a line-imaging VISAR. 
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TABLE I: Parameters of the experiments with Al [cases (a)- 
(d)] and Fe [cases (e)-(h)] targets (see Fig. [I]). Eo is the ion 
energy, No is the number of ions per pulse, r is the pulse 
duration, and Ax x Ay are the FWHM dimensions of the 
beam focal spot. The depth of the beam in the absorber is 
Ar — 0.2 mm in all the cases. All spatial dimensions listed 
are in millimeters. 
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100 
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0.05 
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10 11 


100 


0.3 x 1.0 


0.5 


0.1 
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10 10 


50 


0.3 x 1.0 


0.5 


0.075 



preaching 1 Mbar shall become available with the com- 
pletion of the high-current upgrade of the SIS- 18 heavy 
ion synchrotron of the GSI Hclmholtzzcntrum fur Schw- 
erionenforschung in Darmstadt |31j in the year 2009. In 
RWL experiments with intense ion beams, the typical 
compression times from zero to the maximum pressure 
are about 20 ns for the aluminium and about 10 ns for 
the iron samples considered here, which is comparable 
with the dissipative relaxation times of aluminium and 
iron |32l 133] . Resent laser-driven experiments |23j have 
demonstrated that at such short time scales there exists 
a stiffer response than had been expected from previous 
slower ramp compression experiments and from models 
based on either static or shock- wave experiments. 

The layout of the experiment is shown in Fig. [I] The 
ion beam with an elliptical focal spot and Gaussian trans- 
verse intensity distribution heats an absorber slab along 
its surface. Lead is chosen as the absorber material. 
Since the stopping ranges of energetic 238 TJ ions in lead 
absorber are much larger than the length of the absorber 
slab, the latter is heated uniformly along the z axis. The 
focal spot size and its aspect ratio can be varied, what 
allows for generation of a quasi-planar ramp wave. In the 
similar manner as in the laser-driven ramp compression 
or shock wave experiments [23,, a stepped target with 
rather wide steps along z direction can be used. 

Similarly to laser-driven RWL experiments [2D] , the 
process of the ion beam driven RWL consists of three 
distinct phases. First, the ion beam evaporates the ab- 
sorber material off its surface or unloads the material due 
to a shock wave generated inside the absorber, depending 
on the beam intensity, focal spot size and the distance Ar 
between the beam's axis and the absorber surface. At the 
next stage, the absorber's material expands into the vac- 
uum gap. The adjustable displacement of the beam axis 
relative to the absorber surface allows for some degree 
of control over the plasma expansion regime. In the last 



stage, the absorber's plasma piles up against the target 
foil thus producing a smoothly increasing pressure load 
in the sample. The evolution of the compression wave 
launched into the sample is shown in Fig. [2] One can 
see that the compression wave gradually steepens with 
the distance traveled before it breaks into a shock wave. 
This distance determines the maximal allowed thickness 
of the target foil. 

Provided the compression wave launched into the sam- 
ple does not break into a shock wave, the pressure and 
density history at the front surface of the sample can 
be reconstructed using the bootstrap back integration 
method [13] • The method uses the rear surface veloc- 
ity history as an initial condition for a Riemann solver. 

The major benefit of the proposed experimental design 
as compared to other approaches for EOS measurement 
with heavy ion beams [M] [35] is twofold. Firstly, the 
EOS along a compression adiabat can be determined by 
measuring only one parameter — the rear surface ve- 
locity employing a line-imaging VISAR as the princi- 
pal diagnostics. Secondly, this design does not rely on 
detailed knowledge of the beam-matter interaction pro- 
cesses like the stopping power, since the beam energy is 
not deposited into the sample directly 36J but is con- 
verted to the kinetic energy of the absorber material. 
In order to interpret and use the results of the experi- 
ment one therefore does not have to precisely measure 
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FIG. 2: Evolution of pressure and density for the case (b) 
of Table [I] The dotted vertical lines indicate the location of 
the target front surface at different times. The phases of 
the process are: (1) t = 59 ns — the evaporated absorber 
material — lead is accumulating at the front surface of the 
aluminium sample; (2) t = 78 ns — the isentropic compression 
has started; (3) t — 86 ns — the accumulation continues; the 
compression wave reaches the rear surface of the sample before 
the shock wave has formed; (4) t — 94 ns — the motion of the 
rear surface of the target foil can be detected. 
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FIG. 3: Rear surface velocity and front surface pressure his- 
tories for different cases described in Table |U 



the transverse distribution of the beam intensity at the 
focal plane, which can be problematic for intense focused 
heavy ion beams [37]. Moreover, the presented design has 
also an important advantage that the target is not being 
preheated by energetic secondary particles and projectile 
fragments produced during the interaction of the beam 
with the absorber material. Furthermore, this approach 
requires neither high accuracy of the beam-target align- 
ment, nor good short-to-shot reproducibility of the beam 
parameters. 

Parameters of the beam and the target geometries con- 
sidered in this work are summarised in Table [I] The cases 
(a),(e) and (b),(f) in the table correspond to the beam 
parameters of the SIS- 100 synchrotron to be available 
at FAIR [3B], the cases (c),(g) correspond to the beams 
that will be available after completion of the SIS-18 up- 
grade |31J . and the cases (d),(h) correspond to the present 
SIS-18 beam available for high energy density physics 
experiments [37] • A lead absorber slab of 400 /im 
thickness was considered for all the analysed cases. The 
geometric parameters of the target have been adjusted 
to guarantee the planarity of the ramp compression wave 
propagating in the x direction for at least y < ±50 /im, 
and to ensure a shockless compression of the target. The 
planarity of the designed experiment justifies a one di- 
mensional analysis using the ID hydrodynamic simula- 
tion code [3D] . 

The calculated histories of the pressure at the front and 
the velocity at the rear surfaces of the aluminium sam- 
ples are shown in Fig. [3] The planarity of the compression 
wave can be adjusted by varying the ion beam extension 
in the y direction. However, this also affects the level of 
specific energy density deposited by the beam in the ab- 
sorber. Two possible beam cross sections are compared: 



in the case (a) the FWHM height of the beam spot is 
0.5 mm, whereas in the case (b) — 1mm. Increasing the 
height of the beam spot allows one to increase the vac- 
uum gap which enhances the smoothness of the loading 
without loosing in the planarity. The curves correspond- 
ing to these two cases indicate that the enhancement of 
the smoothness of the loading compromises the ampli- 
tude of the compression: the maximum pressure corre- 
sponding to case (a) is about 3.5Mbar, whereas in case 
(b) it is only about 2 Mbar. The maximum pressure that 
one can obtain using upgraded SIS-18 beam (c) is about 
IMbar and the currently available SIS-18 beam (d) can 
isentropically load aluminium up to 500 kbar. 

The strain rate is approximately 8 x 10 7 s _1 at the 
peak pressure of about 2 Mbar in the examined case (b) 
(see Table [T]) . The corresponding time to the maximum 
compression is approximately 20 ns, which is compara- 
ble to the rise time of low-stress steady shock [32 . The 
laser RWL experiments by Smith et al. 23J at similar 
compression times and strain rates have demonstrated a 
stiffening of the stress-strain response. Using ion beams, 
the adjustable ion pulse duration, focal spot size and the 
depth of the ion beam in the absorber allow a degree of 
control over the compression time, amplitude and strain 
rate. This provides a unique tool to carry out parametric 
studies of the stress-strain response. 

A similar analysis has been carried out for iron sam- 
ples (Table [I]) . The corresponding histories of the front 
face pressure and rear face velocity are shown in Fig. [4] 
One can see that a Mbar pressure loading is within the 
rich of the existing SIS-18 beam, whereas the upgrade of 
the machine can provide the compression data of up to 
about 3 Mbar in the nearest future. Up to 6 Mbar pres- 
sures will be generated in the corresponding experiments 
at FAIR with extremely high strain rates of the order of 
Hps" 1 . 

The curves (e) and (f ) in Fig. [4] demonstrate the dif- 
ferences between the velocity traces in case of shock and 
adiabatic compression, respectively. In order to ensure 
shockless loading, the sample's thickness had to be re- 
duced from 100 /jm as in the case (e) to 50 /im as in the 
case (f). Obviously, decreasing the sample's thickness 
causes the motion of the sample's rear surface to start 
earlier. This motion acts to release the pressure of the 
absorber's material piling up at the front surface of the 
sample and to reduce the pressure amplitude in the tar- 
get. Therefore, a compromise should be devised to ensure 
the adiabatic character of the compression on one side, 
and to obtain the highest possible pressure on the other. 

The capabilities of the high-energy heavy ion beam 
accelerators available at GSI and later at FAIR to drive 
planar isentropic compression of solids were analysed. It 
was shown that the time scales of the loading are com- 
parable with the fastest laser drives ^ZU\ [53], and the 
amplitudes of pressure surpass those obtained using the 
high-power magnetic drives |17j . These features allow to 
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FIG. 4: Rear surface velocity and front surface pressure his- 
tories for different cases of Fe samples described in Table [I] 



investigate the dynamical response of solids in the regime 
of previously unattainable parameters and possibly the 
dynamics of the pressure-induced structural phase trans- 
formations [25] [29, 30J . Moreover, ramp loading with 
pressures higher than 0.1 Mbar can be obtained in multy- 
layer targets |30j using the currently existing SIS-18 ring 
for the purpose of studying the pressure induced phase 
transformations in water. 

Support from EPSRC and INTAS is gratefully ac- 
knowledged. 
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